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Abstract 

We develop a fuhy covariant, weU-posed 5D effective action for the 6D cascading gravity brane- 
world model, and use this to study cosmological solutions. We obtain this effective action through 
the 6D decoupling limit, in which an additional scalar degree mode, tt, called the brane-bending 
mode, determines the bulk-brane gravitational interaction. The 5D action obtained this way 
inherits from the sixth dimension an extra tt self-interaction kinetic term. We compute appropriate 
boundary terms, to supplement the 5D action, and hence derive fully covariant junction conditions 
and the 5D Einstein field equations. Using these, we derive the cosmological evolution induced on 
a 3-brane moving in a static bulk. We study the strong- and weak-coupling regimes analytically 
in this static ansatz, and perform a complete numerical analysis of our solution. Although the 
cascading model can generate an accelerating solution in which the tt field comes to dominate at 
late times, the presence of a critical singularity prevents the tt field from dominating entirely. Our 
results open up the interesting possibility that a more general treatment of degravitation in a time- 
dependent bulk, or taking into account finite brane-thickness effects, may lead to an accelerating 
universe without a cosmological constant. 



I. INTRODUCTION 



A fundamental conundrum exists as to whether the accelerated expansion of the universe 
is due to a new form of energy or novel gravitational physics revealing itself at ultra-large 
scales, extremely low spatial curvatures, and low cosmological densities. Along with studies 
of different forms of dark energy and modifications to gravity, considerable attention has been 
paid to the possible role played by higher-dimensional theories, in which our four-dimensional 
world is considered to be a surface (a "brane" ) embedded in a higher-dimensional spacetime 
(the "bulk"). In the old Kaluza-Klein picture it was necessary for the extra dimensions 
to be sufficiently compact (for reviews see, e.g. [I1I2])- Recent developments, however, are 
based on the idea that all standard model particles are confined to a 4D brane, whereas 
gravity is free to explore the bulk [3HS]. As such, "large" extra dimensions are conceivable, 
giving rise to a much smaller fundamental Planck mass than the effective Planck scale we 
observe today [H El [7] . A well-studied example of such a theory is the Dvali-Gabadadze- 
Porrati (DGP) model [8], in which our observed AD universe is embedded in an infinite fifth 
dimension. In this picture, the higher-dimensional nature of gravity affects the AD brane 
through deviations from general relativity on horizon scales, r ~ cH^^ (where c is the speed 
of light and Hq is the Hubble constant), that may give rise to the observed accelerated 
expansion. 

In the DGP model, integrating out the bulk degrees of freedom yields an effective action 
for the AD fields containing, besides the graviton, an extra scalar degree of freedom, vr, 
called the brane bending mode [914TT]. The vr field contributes to the extrinsic curvature 
of the boundary and interacts strongly at the energy scale A5 = MI/M4, where M5 and 
M4 are the 5D and AD Planck masses respectively. In analogy with massive gravity [T2] . 
there exists a decoupling limit in which the strong interaction scale A5 is held fixed while 
M4, M5 — )■ 00. All other degrees of freedom (including the graviton and a vector A*"^) decouple 
in this limit. This implies that the dynamics of the scalar field vr can completely describe 
all interesting features of the DGP model, including the Vainshtein screening effect [H] 
and the self-accelerated cosmological solution [13]. It has now been established that the 
branch of solutions that include self-acceleration suffers from ghost-like instabilities flU\ \1M - 
[T8] . On the observational front, DGP cosmology is statistically disfavored in comparison 
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to ACDM [T9H2T] and is significantly discordant with constraints on the curvature of the 
universe |22j . 

Recently, a phenomenological approach to the cosmological constant problem — degravita- 
tion [231425] — has been developed. In degravitation it is postulated that the cosmological 
constant is indeed responsible for dark energy. The cosmological constant problem that the 
observed value is at least 120 orders of magnitude smaller than vacuum energy density pre- 
dicted theoretically, is solved not by making the vacuum energy density small, but instead, 
by having a large cosmological constant whose gravitational effect is suppressed by making 
gravity extremely weak on large scales. The DGP model should, in principle, provide a more 
fundamental implementation of degravitation. However, the weakening of gravity observed 
in DGP is insufficient to account for the disparity between the expected and observed values 
of the cosmological constant. This fact, in addition to the above mentioned problems of the 
DGP model, have led to the idea of cascading DGP [26ti30] — a higher-dimensional gen- 
eralization of the DGP idea, which is free of divergent propagators and ghost instabilities. 
In this model one embeds a succession of higher-codimension branes into each other, with 
energy-momentum confined to the AD brane and gravity living in higher-dimensional space. 
(See [3T] for a related framework.) 

The implemention of degravitation within the cascading gravity idea provides an intriguing 
new theoretical avenue for solving the problem of dark energy. However an important litmus 
test is whether such models can reproduce a successful cosmological evolution. Studies thus 
far in this direction have assumed an effective AD cosmology for degravitation by generalizing 
that for DGP [321 [33]. However, to perform a more complete study of cosmology on the 
brane, it is necessary to integrate out the sixth and fifth dimensions to obtain a AD effective 
theory. 

In this paper we start from the action for cascading gravity in QD and obtain an effective 
linearized action in the decoupling limit. This gives rise to an extra brane-bending 
scalar degree of freedom (the vr field) in the bD action. As a proxy for the complete 6Z? 
cascading set-up, we propose a bD non-linear and covariant completion of the quadratic 
action. A similar strategy was used in [31], where an analogous AD covariant action was 
shown to reproduce much of the phenomenology of the full DGP model. In our case, the 
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resulting action is a 5D scalar-tensor theory, describing 5D gravity and a scalar vr, coupled 
to a 4D brane. Because of its scalar-tensor nature, the standard Israel junction conditions 
must be revisited. We derive the appropriate junction conditions across the 4D brane using 
two different techniques. These can then be used in conjunction with the bulk equations 
to study cosmology on the brane. For concreteness, we consider the cosmology induced on 
a moving brane in a static bulk geometry. We find analytical solutions in the strong- and 
weak-coupling regime for the vr field, and numerically integrate the full equations of motion. 
Thanks to the Vainshtein screening mechanism, the resulting AD cosmology is consistent 
with standard big bang expansion history at early times, but deviates from ACDM at late 
times. We find that tt contributes to cosmic acceleration at late times, but a singularity in 
the brane embedding prevents vr from accounting for all of dark energy. 

In sectionjnjwe outline the 6D cascading gravity model we consider, and propose an effective, 
covariant 5D action with a strongly interacting vr field that encodes the 6D physics. In 



section III we derive the appropriate boundary terms necessary in order for our action to 



have a well-defined variational principle. The resulting bulk equations of motion and brane 



junction conditions are computed in section IV We then turn in section M to the search for 



cosmological solutions on the brane, by considering its motion in a static bulk. Finally, we 



draw together our findings and discuss implications in section VI 



A comment on our notation: we denote coordinates in the full six dimensional spacetime by 
X ^ Ob ^ Ob ^ OC ^ OC ^ 00 ^. Indices M,N, ... run over 0,1,2,3,5 (i.e. the 4 + ID coordinates), indices 
/i, z/, ... run over 0,1,2,3 (i.e. the 3 + ID coordinates), and indices ... run over 1, 2, 3 (i.e. 
the 3D spatial coordinates). We further denote the fifth and sixth dimensional coordinates 
hj y = and z = x^, where convenient. 



II. A PROXY THEORY FOR CASCADING GRAVITY 



The DGP model consists of a 3-brane embedded in a fiat, empty 4 -|- 1-dimensional bulk. 
Despite the fact that the extra dimension is infinite in extent, the inverse-square law is 
nevertheless recovered at short distances on the brane due to an intrinsic, four-dimensional 
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Einstein-Hilbert term in the action 

r f /M^ \ 

St,GP = / d^Xv^^i?5 + / d^Xv^ -fRA + /:matter ■ (l) 
^bulk ^ ^ branc V ^ / 

The Newtonian potential on the brane scales as 1/r at short distances, as in 4D gravity, 
and asymptotes to at large distances, characteristic of 5D gravity. The cross-over scale 
rrig ^ between these two behaviors is set by the bulk and brane Planck masses via 

From the point of view of a brane observer, this force law arises from the exchange of a 
continuum of massive gravitons, with setting an effective mass scale for gravity on the 
brane. The DGP model is therefore a close phenomenological cousin of Fierz-Pauli massive 
gravity. In particular, brane gravitons form massive spin-2 representations with 5 helicity 
states, with the helicity-0 mode having a small strong-coupling scale, 

A, = {mlM,Y/'. (3) 



There are many reasons to consider extending this scenario to higher dimensions: 



• Pragmatically, cosmological observations already place stringent constraints on the 
DGP model [T^M^ . In higher dimensions, however, the modifications to the Fried- 
mann equation are expected to be milder, which traces back to the fact that the AD 
graviton mass term is a more slowly-varying function of momentum. The resulting 
cosmology is therefore closer to the ACDM expansion history, thereby allowing a wider 
range of parameters. 

• Another motivation, as we have already mentioned, is the degravitation idea [2H [25] 
for addressing the cosmological constant problem; namely that gravity acts as a high- 
pass filter that suppresses the contribution of vacuum energy to the gravitational 
field. Although the infrared weakening of gravity displayed in the DGP force law is 
suggestive of a high-pass filter, in practice this weakening is too shallow to "filter out" 
vacuum energy. However, the situation is more hopeful in D > 5 dimensions, where 
the force law on the brane falls more steeply as l/r^~^ at large distances [25] . 
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While seemingly a straightforward task, generalizing the DGP scenario to higher dimensions 
has proven challenging. To begin with, the simplest constructions are plagued with ghost 
instabilities around flat space [551 155] . Another technical hurdle is the fact that the 4D 
propagator is divergent and requires careful regularization [371 EH]- Finally, for a static 
bulk, the geometry for codimension N > 2 has a naked singularity at a finite distance away 
from the brane, for an arbitrarily small tension [23] . 

It was recently shown that these pathologies are absent if the 3-brane is embedded in a 
succession of higher-dimensional DGP branes, each with their own Einstein-Hilbert term. 
In the 5 + 1-dimensional case, for instance, the 3-brane lies on a 4-brane, with action. 



As a result, the force law on the 3-brane "cascades" from to to as one moves 
increasingly far from a source, with the AD — )■ 5D and 5D — )• 6D cross-over scales given 
respectively by mg ^ and mg \ with 

This cascading gravity setup is free of the aforementioned pathologies: the theory is pertur- 
batively stable provided that the 3-brane is endowed with a sufficiently large tension [261 EZ] ; 
the 5D Einstein-Hilbert term acts as a regulator for the induced propagator on the 3-brane; 
and, as has been shown explicitly for D = 6,7, adding tension on the 3-brane results in 
a completely smooth bulk geometry (except of course for the delta-function singularities 
at the brane locations) and leaves the 3-brane geometry flat, at least for sufficiently small 
tension [28] . 

The next question is, of course, whether the resulting cosmology is consistent with current 
observations, and, more interestingly, whether it offers distinguishing signatures from ACDM 
cosmology. Unfortunately, finding analytical solutions is a hopeless task, even in the sim- 
plest QD the bulk metric is generally expected to depend on all extra-dimensional 
coordinates plus time [5U] . 

In this paper, we instead study a 5D "proxy" brane-world theory for 6D cascading gravity, 
consisting of a scalar-tensor theory of gravity in the 5D bulk. This is obtained by generalizing 
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the well-known decoupling limit of standard DGP |9] to the cascading case. The limit we 
propose is M5, Mq — > 00, with the strong-coupling scale 

kept fixed. In this limit, the action (|4]) may be expanded around fiat space, and reduces to 
a local theory on the 4-brane, describing 5D weak-field metric perturbations Hmn and an 
interacting scalar field tt. The latter is the helicity-0 mode of massive gravity on the 4-brane, 
and has a geometrical interpretation as measuring the extrinsic curvature of the 4-brane in 
the 6D spacetime. The resulting action is [26] 



^decouple 2 "-l 



bulk 

+ / dV 

' brane 



L, J.D77Zg 

1 



(7) 



where 



{EK)mn — —^{O^hMN — VMN^^h — OmO^ hxN 

- dNd^huK + VMNd^d^hKL + dudNh) (8) 

is the linearized Einstein tensor in SD, and {Eh)^^ that in AD. To see that only these terms 

Q In 

survive in the decoupling limit, introduce canonically-normalized variables tt'^ = M5 vr and 
^MN — ^5 hMNi which have the correct mass dimension for scalar fields in 4-1-1 dimensions. 
The quadratic terms in ([T]) become independent of M5 under this field redefinition, whereas 
the cubic term reduces to {dTXc)'^^^TXc/ ^J"^ ■ All other interactions in (ij are suppressed by 
powers of I/M5, l/Mg and therefore drop out in the decoupling limit. 

In using ([T]) as our starting point, we are motivated by the fact that nearly all of the 
interesting features of DGP gravity are due to the helicity-0 mode vr and can be understood 
at the level of the decoupling theory [101 El]. Of course, as it stands ([T]) is restricted to 
weak-field gravity and is therefore of limited use for cosmological solutions. As our "proxy" 
brane- world scenario, we propose to complete ([T]) into a covariant, non-linear theory of 
gravity in bD coupled to a 3-brane. By construction, the weak-field limit of our theory will 
coincide with ([t]). A similar approach was followed in [M] to mimic the hD DGP scenario 
with a proxy effective theory in 4D. Despite being a local theory in 3 -|- 1 dimensions, the 
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resulting cosmology was found to be remarkably similar to that of the full 4 + 1-dimensional 
DGP framework, both in its expansion history and evolution of density perturbations. 

Generalizing the strategy of |3lj to the cascading gravity framework, we are led to propose 
the following non-linear completion of ([T]): 



bulk 



J brane 



2 -Tl^A ~r *— matter 



(9) 



It is straightforward to check that this theory indeed reduces to ([t]) in the weak-field limit, 
and therefore agrees with cascading gravity to leading order in I/M5. (This is most easily 
seen by working again with the rescaled variables vr^ and h\j^.) The proposed completion 
is by no means unique, since one could consider a host of Ms-suppressed operators which 
would disappear in the weak-field limit. Our hope is that the salient features of cascading 
cosmology are captured by our effective theory, and that the resulting predictions are at 
least qualitatively robust to generalizations of (|9|. 

The effective action ([9]) must be supplemented with suitable boundary terms in order to 
yield a well-defined variational principle. Other than a Gibbons-Hawking- York-like term, 
the form of the cubic term in vr clearly necessitates its own boundary contribution. In the 
next section, we derive these boundary terms, which will be essential in deriving the junction 
conditions. 



III. BOUNDARY TERMS IN THE bD EFFECTIVE THEORY 

Because of the form of the cubic term, varying (|9]) with respect to tt yields contributions 
on the 3-brane of the form ~ (97r)^£„(57r, where £„ is the Lie derivative with respect to the 
normal. Such terms cannot be set to zero by the usual Dirichlet boundary condition, ^vr = 0, 
and must therefore be canceled by appropriate boundary terms in order that the action be 
truly stationary and the variational principle be well-posed. Gravity also requires its own 
boundary contribution, which is a generalization of the well-known Gibbons-Hawking- York 
term |10l HI]. (We should, of course, also include boundary terms at infinity, but we will 
ignore these since they do not play any role in the junction conditions.) 
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To derive the boundary terms, it is convenient to work in the Arnowitt, Deser and Misner 
(ADM) coordinates [42j, with y playing the role of a "time" variable, 

ds%) = N^dy^ + g^,(dx^ + mdy){dx'' + N'^dy) , (10) 

where N and A^^ are the lapse function and the shift vector, respectively. In the "half- 
picture", the bulk extends from ?/ = to oo, and the 3-brane is located at ?/ = 0, with 
normal vector um = (0, 0, 0, 0, A^). 

In ADM coordinates, the 5D Einstein-Hilbert term takes the form 

^gravity = ^ / d^xdy v^ATg-^-Z^ + _ j^^^j^.u ^ 2V,, (n^V^n*^ - n^'K)] , 
-'y>o 

(11) 

where K^^^, is the extrinsic curvature tensor 

K^, = ^£„g^, = ^ {dyq^, - D^N, - D,N^) . (12) 
Here is the covariant derivative with respect to the 4Z) induced metric q^y. Unlike 



^-gravity = ^ / d^xdy v^A^ [e'^^/^ (i?4 + ' K^^K^'' - 3KCnn) - 2046-3^/2] 



standard gravity, the CnK term in (11) is not a total derivative and must be treated with 
care. Integrating by parts gives 

+ Ml ! d^x^qt-^^l'^K , (13) 
and the last term must therefore be canceled with a Gibbons- Hawking- York boundary term 

A^GHY = -Ml [ d^x^qe-^^l'^K . (14) 

Jj/=0+ 

Similar considerations for the 7r-sector lead us to require adding the boundary term 

27 /■ / 1 \ 

= -^^ d^xv/^ d.nd^nCnn + - (£„7r)3 , (15) 

6Z rriQ Jy=Q+ \ 6 / 

where 

CnTT = N-\dy - N^'d^)!: . (16) 

Note that in the flat space limit this agrees with the vr boundary term obtained in in 
the decoupled theory. 
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Including (14) and (15), the full 5D action is therefore 

27 



S 



Ml 



bulk 



Ml 



+ 



brane 
4 



d Xa/^ 



brane 



2 



IQml 
27 



32ml 



Ra ~\~ -^matter 



(17) 



Although we obtained the boundary terms using the ADM formalism, the result is fully 
covariant and hence holds in any coordinate system. In particular, given the unit normal 
vector to the brane n*^ in a general coordinate system, the Lie derivative is given by ^nvr = 
n'^'^dM'^, and the induced metric by qmn = Qmn — ^m^n- One can check that varying this 
action with respect to the metric and vr does not yield any normal derivative terms of the 
form CnSq^i, and CnSir on the boundary. 



IV. COVARIANT EQUATIONS OF MOTION ON AND OFF THE BRANE 



Our goal now is to derive the bulk equations of motion and brane junction conditions that 



result from (17). (See [UHlTj for earlier work on junction conditions in scalar-tensor brane- 



world scenarios.) Starting with the bulk, varying (17) with respect to the metric yields the 
Einstein equations 



27 



MN 



9(Af(97r) dN)TT - -QMNd (dir) Oktt - dAfirdNirn^ir 



16m| 

{gMND5 - Va/Vjv) e-3"/2 



(18) 



where Gmn the 5D Einstein tensor. The second line is typical of scalar-tensor theories and 
arises from the non-minimal coupling of vr to gravity. Varying with respect to vr, meanwhile, 
gives 

4 



(19) 



where Rmn is the Ricci tensor and -R5 is the Ricci scalar. Remarkably, even though 



the cubic vr interaction in (17) has four derivatives, all higher-derivative terms cancel in 



the variation, yielding a second-order equation of motion for vr. This is a nontrivial and 
important property of the DGP tt lagrangian [10]. In the decoupling limit of Fierz-Pauli 
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massive gravity, by contrast, the vr lagrangian takes an analogous form, but its equation of 
motion is higher order — there is a ghost mode propagating at the non-hnear level 1^51 - 
\5U\ . See [SH [S2] for an interesting recent proposal of a non-linear completion of Fierz-Pauli 
gravity that seemingly avoids these pitfalls. 

Next we obtain the junction conditions at the brane position by setting to zero the boundary 



contributions to the variation of (17). Assuming a Z2-symmetry, variation with respect to 



the metric yields the Israel junction condition 



27 Ml 



» me V 3 J 

+ Tl^J - MlG^^ , (20) 



where 



y(4) 



2 5{^/^C 



matter ) 



(21) 



is the matter stress-energy tensor on the brane, and is the Einstein tensor derived from 



the induced metric q^y. Similarly, varying (17) with respect to the scalar, we obtain after 
some algebra the boundary condition for vr on the brane: 
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OTTla, V / 



(22) 



Equations (20) and (22) are not independent, of course; the divergence of (20) can be 



shown to be proportional to (22) after using the bulk momentum constraint equation. As 



a nontrivial check on our junction conditions, we have evaluated (20), (22) in a gauge in 



which the brane is at fixed position {y = 0) and the bulk metric is time- dependent, and have 
shown that the result agrees with the boundary conditions obtained by integrating the bulk 



equations (18)-(19) across the delta-function sources at y = (see Appendix 



V. THE COSMOLOGICAL EVOLUTION ON THE BRANE 

The study of brane-world cosmology requires us to use our equations of motion to obtain a 
Friedmann equation on the brane, assuming homogeneity and isotropy along the 3+1 world- 



volume dimensions. The junction conditions (20) and (22) do not form a closed system of 



equations for q^^, hence deriving an induced Friedmann equation requires knowledge of the 
bulk geometry [39]. 
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Because of the bulk scalar field, there is no Birkhoff's theorem to ensure that the bulk 
solutions are necessarily static under the assumption of homogeneity and isotropy on the 
brane — the most general bulk geometry depends on both the extra-dimensional coordinate 
and time. For concreteness, however, we focus here on a static warped geometry with 
Poincare-invariant slices, 

dsl,^ = a\y)i~dT' + dx') + dy\ (23) 

While admittedly restrictive, we view this ansatz as a tractable first step in exploring cas- 
cading cosmology. And, as we will see, the resulting phenomenology is already surprisingly 
rich. 

The brane motion is governed by two functions, y{t) and T(t), describing the embedding, 
where t is proper time on the brane. The induced metric is of the Friedmann-Robertson- 
Walker (FRW) form, with spatially-fiat {k = 0) constant-time hyp ersurf aces, 

dsLne = -d^' + a'(y)df^ (24) 

where, by virtue of t being the proper time. 



Given a solution a{y) to the bulk equations ([18|)-([19|), the covariant junction conditions (20) 



and (22) allow us to solve for the embedding (?/(t), r(t)), and hence the cosmology induced 



by brane motion through the warped bulk. 



A. A Dynamic brane in a static background 



With the static ansatz (23), the bulk equations (18)-(19) take on a form reminiscent of 
cosmological equations, with a{y) acting as a scale factor as a function of "time" y. In 
particular, the (5, 5) component yields a Friedmann-like equation 
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V a , 





1 



whereas the (yU, v) components yield 

,/\ 2 



a 
a 



+ 



9 



16mg 



e ' TT TT 



-TT 



/2 



TT 



(26) 



(27) 
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Meanwhile, the equation of motion for vr can be written as 

2 A 



d 

dy 



a 



-7r'2 1 +4 



a'n' 



--rule-'-!- 



(28) 



As usual, the Bianchi identity guarantees that only two of these equations are independent. 
Finding exact solutions to these equations requires a numerical approach, which we will 



perform in section VP To offer analytical guidance, however, we seek approximate solutions 
to (26)-(28) in the so-called strong- (section [V B ) and weak-coupling (section VC) regimes 
in which the non-linear terms in vr respectively dominate or are negligible in these equations. 

The brane embedding (?/(t),r(t)) is determined by the junction conditions, which involve 



the extrinsic curvature tensor and the Lie derivative of vr. Using (25) the relevant quantities 
are 

' " 1 d 

a dy 



1 + 



d^ 
dt 



3 ' 



if 




(29) 



and 



£n7r 



TT 



1 + 



dy 
dt 



(30) 



For the stress energy on the brane, we assume a collection of (non-interacting) perfect fluids 
with energy densities pm and pressures .Pm , obeying the standard continuity equations 

dpm 



dt 



3i/(p» + P«) = 



(31) 



where if = dlna/dt is the Hubble parameter on the brane. These components may include 



baryonic matter, dark matter, radiation and a cosmological constant A. Equation (31) is 



consistent with the picture that matter is not allowed to flow into the bulk and is confined 
to the brane. 

It is clear, therefore, that given a bulk solution a{y), a single junction condition is sufficient to 



solve for the cosmological evolution on the brane. Indeed, although (20) and (22) yield three 



equations, two of these follow from the bulk Hamiltonian and momentum constraints, which 
are automically satisfied given a solution a{y). Since we are interested in the Friedmann 



equation on the brane, the natural choice is the (0,0) component of (20). Noting that 
(9o7r = Tx'dy/dt and dy/dt = aH/a', we can write the resulting equation as the standard 
Friedmann equation with an additional effective energy density resulting from the tt 
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field, 

?,H^Ml = ^ p« + , (32) 

i 

where 

p. . M|V;OT{3i,(2(^)^I)^-e.-(|:4)} , (33) 

encoding all the complexity and new physics of our model. Given a solution a{y), 7r(y) 
to the bulk equations, we may invert this relation to obtain y{a), and use this to express 



all ?/-dependent terms in as functions of a. Equation (32), together with the continuity 



equations (31), then form a closed system for the brane scale factor a{t). 



Before moving on to explicit solutions, we note in passing that p^^ is not positive definite. 
When combined with A, this can lead to an effective equation of state parameter w < —1 for 
the effective dark energy component. This phantom behavior already occurs in the normal 
branch of the standard DGP model [531155] . a phenomenon that can be understood in the 
decoupling limit as arising from non-minimal coupling of the brane-bending mode to brane 
gravity [33]. (It is well-known that w < —1 can be achieved in scalar-tensor theories when 
working in the Jordan frame [561 - I59] .) Similarly here, the scalar vr is kinetically-mixed with 
the brane graviton, which can lead to phantom behavior for dark energy. 



B. The Strong-coupling regime 



By analogy with the Vainshtein screening mechanism around astrophysical sources, we ex- 
pect that at early times, when the energy density in the universe is high, vr should be strongly 
coupled and cause small deviations from standard AD Friedmann cosmology. In other words, 
the non-linear terms in vr dominate, but as a result is negligible compared to matter and 
radiation. Moreover, since the total variation in vr is expected to be small in this regime 
(|A7r| ^ 1), by rescaling M5 we can assume that e^^/"^ ~ 1. 
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Consider (26) and (27) in the regime in which the non-hnear terms in vr dominate: 

a' 9 



TT 



a" fa'V 9 



+ [-) = TT^vr'^vr" . (34) 
a \ a / Ibmi 
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These admit scahng solutions, given by 

a{y) = 

^(y) = \^\li^e\y\ j , (35) 

where we have a chosen a mass scale proportional to rriQ in the solution for a{y). This leaves 
the scale factor today, oq, to be a free parameter. It is straightforward to check that the above 




solution also satisfies the third bulk equation (28 ) in the strong-coupling approximation. The 
approximation vr ^ 1 implicit in (35) is therefore valid provided y ^ mg^. This defines the 



regime of validity of this solution. 

The naked singularity at y = — the analogue of a big bang singularity in cosmology — 
introduces a plethora of complications if included as part of the bulk geometry. It is therefore 
safest to exclude this part of the geometry when performing the Z2 identification. As a result, 
however, the warp factor grows without bound as one moves away from the brane, which 
may indicate a strong-coupling problem. A related question concerns the stability of this 
solution — by analogy, the self-accelerated branch of the DGP model also has a growing 
warp factor [T3] and is well-known to suffer from instabilities. We leave a careful study of 
these important issues to future work. 

The above solutions for 7i{y) and a{y) can be used to express the effective Friedmann equa- 
tion (32) solely in terms of the brane scale factor. In the strong-coupling regime, the vr'^/mg 



term dominates over the vr' term in (33), giving 



,,^Ml{2^^a^'^^,]f^7{^\ (36) 



and (32) reduces to 
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where is defined in ([2]). Combined with the matter fluid equation (31), this effective 
Friedmann equation completely describes the evolution of the universe in the strong-coupling 
regime. 

In contrast with the standard DGP Friedmann equation, H'^ = p/3M| =p 2m^H , where 
the departure from AD gravity is set by H/m^, here the relative importance of p^^ also 
depends on a time-dependent scale rriQ/a}'^/^. In particular, for a flxed initial value of a, 
the magnitude of the modiflcation can be set arbitrarily by a suitable choice of itlq. This 
freedom reflects the choice of initial condition for the brane motion in the bulk — because 
the bulk is warped, different initial locations of the brane yield different expansion histories. 
In the standard DGP model, on the other hand, the bulk is flat Minkowski space, and hence 
all initial conditions (within the same branch of solutions) are related by the Poincare group. 

To proceed, we consider two limiting cases: 



\i H ^ niQ/a}'^^^ ^ then the modiflcation to the Friedmann equation further reduces to 



^ 2Ml — . (38) 

Assuming that the universe is dominated by a matter component with general equation 
of state w, then H ~ a^^^^^'^^^'^, and thus ~ a3(i-i5«))/io^ terms of an effective 
equation of state for the tt fleld, deflned through dlnp^r/dlna = —3(1 -|- Wt,), we have 

11 3 

^- = -^^ + 2^- (39) 

In particular, since w-j^ < w, it is clear that p^^ becomes more and more negligible as 
we look backwards in time. Moreover, in a universe dominated by baryonic and/or 
cold dark matter {w = 0), the n fleld can act as a dark energy fluid with phantom 
equation of state w-„ = —11/10. 

A phantom equation of state opens up the possibility of the vr fleld acting like dark 
energy and driving cosmic expansion. In the strong regime, the Friedmann equation 



(32) can be approximated by the cubic equation. 



AH' -3H^ + -^ = (40) 
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with p^/M| = AH^, and A = 2a?'^/^m^/m1. Differentiating (40) gives 

uw--l-3j^ = -l-3^^. (41) 

For fi^ > 5/24 ^ 0.21, this gives w^q < —1/3, and acceleration occurs. However, the vr 
field is unable to dominate the energy density and fully account for the current phase 
of accelerated expansion, because of a singularity at fi^ = 2/3 for which — — oo. 

In the opposite regime, H <C m^^ja}^^!^ ^ we have 

/^^-SMl-^. (42) 

In this case, the vr component has a fixed effective equation of state, = —1/5, 
independent of the matter on the brane. Again, this pushes the total equation of state 
to more negative values. 



C. The Weak-coupling regime 



By analogy once again with the Vainshtein story in DGP, at late times we expect the non- 
linear terms in vr to be negligible, corresponding to gravity becoming higher-dimensional. In 



this approximation, the bulk equations (26) and (27) reduce to 



a 
a 



IT 



+ 



2 V 2 



'2 



a'n' 



-TT ■ 



TT 



(43) 



which again admit a scaling solution 



12 



2/5 



\n{m6\y\) . 



(44) 



The mass scale in the solution for a{y) has been chosen to be consistent with the strong- 
coupling solution. It is straightforward to check that this solution consistently satisfies the 



third bulk equation (28) in the weak-coupling limit. 
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Substituting this solution into (33), the effective energy density in vr in the weak-couphng 
regime reduces to 



In the hmiting case in which H ^ me/a^/^, this further reduces to 

9 Ml 



p^^^H', (46) 
which imphes that 

^n = ^ + ^w. (47) 

It is interesting to note that for a cosmological constant with w = —1, the vr field also 
behaves as a cosmological constant, = —1. Similarly, for H <^ rrie/a^^'^, 

which behaves like a relativistic component {w.„ = 1/3) independent of the matter on the 
brane. 



D. Numerical Solutions 



To complement the analytical strong and weak-field limits in V B and V C , we numerically 



evolve the full bulk and brane equations given in (26)-(28) and (31)-(33), in the presence of 



matter on the brane. We assume zero spatial curvature on the brane, and include relativistic 
and pressureless components consistent with the standard cosmological model: Q^a = 0.3, 
Qr = 8.5 X 10""^. We further fix the scale factor today to be oq = 1. 

Starting well into the radiation dominated era, with a <C 1, we evolve tt, tt', y and t forward 



with respect to In a: (27) and (28) combine to form an equation for vr", from which we form 



an equation for dvr'/dlna = it" /{a' /a), and (32) can be rewritten as a cubic equation in H^, 
which, if a positive, real solution exists, can be used to evolve t through dt/dlna = 1/H. 

In Fig. [T| we show numerical confirmation of the analytical dynamical attractor solutions 



for discussed in VB and V C For scenarios with tjiq ^ if we find attractor solutions 
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a 



FIG. 1: Evolution of the effective equation of state, = —1 — (l/3)dlnp7r/dlna, for the tt- 
dependent modifications to the Priedmann equation ( |33[ ). The numerical results are consistent 
with the analytical predictions for the large (upper panel) and small (lower panel) mg limits in the 
strong- (a <^ 1) and weak-coupling (a ^ 1) regimes. Here we use the numerical values (in natural 
units c = h=l): Hq = 2.33 x 10^^ Mpc~^ (i.e. Hq = 70 km s^^Mpc'i), (upper panel) = 10^° 
Mpc~^ (mg ^ H) and = 10"^'^ Mpc~^ and (lower panel) mg = 10^^^ Mpc^^ (mg <^ H) and 
ms = 10^^^ Mpc^^. The vr field is a subdominant component of the total energy density at all 
times, and late-time acceleration is driven by A. 

of Wt, = —0.6 and —1.1 in the (strongly-coupled) radiation and matter dominated eras 
respectively, and —1 in the (weakly-coupled) A-dominated epoch. For mg ^ H, strongly- 
and weakly-coupled attractors arise with = —0.2 and = 1/3, respectively. 

The TT field has an effective 'phantom' equation of state in the matter-dominated (mg ^ 
H) regime. This opens up the apparent possibility of cosmic acceleration arising within 
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cascading cosmology without the need for a cosmological constant. However, as discussed in 



Sec. VB , while it is possible to generate acceleration at late times, one hits a singularity in 
the expansion history when = STrGpT^/SH"^ = 2/3 so that the universe cannot smoothly 
transition towards Q.^ 1. In Fig. |2] we show a realization of such a scenario, with the 
onset of cosmic acceleration, and the limiting presence of the singularity. 

This singularity is of an unusual nature — it is not equivalent to the Big Rip scenarios in 
which H and a both become infinite in a finite space of time, since the Hubble parameter H 
and scale factor a remain finite while H diverges. Moreover, the bulk geometry is smooth at 
that point, and it is the brane embedding that is singular. It is possible that this singularity 
could be circumvented by the use of a more general metric ansatz than the static case 
considered here to obtain solutions on the brane, or by accounting for finite brane-thickness 
effects. We leave this to future investigation. 



VI. CONCLUSIONS 



Cascading gravity is a phenomenologically rich framework for exploring new phenomena 
associated with infrared-modified gravity, and offers a promising avenue for realizing de- 
gravitation. This construction circumvents many of the technical hurdles of earlier attempts 
at higher-dimensional extensions of DGP: the induced propagator is free of divergences, the 
theory is perturbatively ghost-free, and adding a small tension on the AD brane yields a 
bulk solution which is nowhere singular and remains perturbative everywhere. Due to its 
higher- dimensional nature, however, extracting cosmological predictions presents a daunting 
challenge. 

In this paper we have considered the more tractable problem of a 5D effective brane-world 
set-up, obtained from the full 6D cascading theory through the decoupling limit. Strictly 
speaking, the decoupling limit leaves us with an action describing a scalar vr and weak-field 
gravity, which is therefore of limited use for studying cosmology. But since vr is responsible 
for most of the interesting phenomenology of cascading gravity, we have proposed a fully 
covariant, non-linear 5D completion of the decoupling theory, as a proxy for the complete 
6D model. Our effective action describes 5D DGP gravity with a bulk vr scalar field, coupled 
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to a AD brane with intrinsic gravity. 

Upon supplementing the 5D action with boundary terms (to yield a well-posed action prin- 
ciple), we obtained covariant junction conditions across the brane, relating the extrinsic 
curvature to delta-function sources on the brane. In order to study cosmology on the brane, 
we then considered a scenario in which a dynamic brane moves across a static bulk, and 
consistently solved the bulk and brane equations of motion. We derived analytical solu- 
tions for the induced cosmology in the strong- and weak-coupling regimes, valid at early- 
and late-times, respectively, and confirmed these expectations with a complete numerical 
analysis. 

Thanks to a cosmological Vainshtein mechanism, the bulk scalar vr and the helicity-0 mode 
of the 4D massive graviton both decouple at early times, resulting in an early-universe 
cosmology that closely reproduces the expansion history of the standard big bang theory. 
At late times, however, these scalar modes effectively contribute to dark energy through a 
modification of the Friedmann equation and result in small deviations from ACDM expansion 
at late times. Although these scalars thus affect dark energy, a singularity in the brane 
embedding prevents the modification from being entirely responsible for cosmic acceleration. 

We are currently studying the evolution of cosmological perturbations in this context. Such 
an analysis should also shed light on the all-important question of stability. With our branch 
choice, the modification to the Friedmann equation behaves as an effective component with 
positive energy density. At first sight this is worrisome, since the counterpart in standard 
DGP is the self-accelerated branch, which is plagued with ghost instabilities. It is crucial to 
investigate whether or not this is the case here too. From a phenomenological perspective, 
we are performing a full likelihood analysis for the predictions of the model, including both 
expansion and growth histories. On much smaller scales, we are also working to derive the 
consequences for Lunar Laser Ranging observations, thereby generalizing the analysis of [60] 
for DGP to the degravitation/cascading framework |25j . 
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Appendix A: Junction conditions for a specific choice of metric 

As a non-trivial check on our covariant junction conditions, in this Appendix we show that 
they agree with direct integration of the bulk equations of motion, when specialized to a 



gauge in which the brane is at fixed coordinate position. Varying the action ( 17) with respect 



to the metric and vr led to the 5D Einstein field equations ( 18 ) and the vr equation of motion 



(19). Upon adding in contributions from delta-function sources on the brane to the bulk 



Einstein field equations we obtain 
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MN 



16mg |_ 



where Tj^^ and G)^J are defined in section 



IV 



The TT equation is as before, 
4 



9 



(A2) 



We are interested in studying brane-world cosmological solutions, for which we specialize to 
5D spacetime metrics of the form. 



(A3) 



with TT = -nij^y). The brane is defined by the hypersurface y = 0, and we compute the 
junction conditions, relating the jump across the brane of a', n\ and vr' or the extrinsic 
curvature to delta-function sources on the brane, for this metric. 
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The metric is required to be continuous across the brane in order to have a well-defined 
geometry. However, its derivatives with respect to y may be discontinuous across y = 0, 
and therefore the second derivatives with respect to y will contain a Dirac delta function 
|61j . so that a" = a" + [a']5{y), and similarly for n. Here a" is the standard derivative 
(the non-distributional part of the double derivative of a), and [a'] is the jump in the first 
derivative across y = Q. We similarly allow vr to be discontinuous across the brane and 
write, vr" = vr" + [7r']5(?/). We further impose a Z2-symmetry {y o —y) across the brane, 
and therefore [X'] = 2X'(0+) for the metric and vr discontinuities. 

In order to obtain the junction conditions we integrate the (0,0) and («, j) components of 



the bulk field equations (Al) and the tt equation of motion ([A2j) over an infinitesimal region 
of the extra dimension y, spanning ?/ = 0. This picks out coefficients of 5{ii) and leads to 
the following three junction conditions, 

a'r, lii'A 9 1 /IttA^ 7r^\ t:' 1 x-^ ^•^ 1 



g-37ro/2 ( "-0 ___(_!10__^^^ 



2 60/ lQml\?,hl nlJ bo GMf^'^"' 2m5 agng 



V + — ^ ( A4) 



riobo 2 bo J 16ml\3bl nljbo 



L f 3 V P« + 2 V o« 1 — (^ + 2^ - 2^^ 



e-o/2(3<. + !^^+9i_,J_![o^ 



(A5) 



ao no J 4m^ |^ 

9 1 r n'n / vrf? TT, 



TTo / ho 


-3 


do 


nl \no 


ao 




<^ 


)• 









+ i^p 7J-^)+3^^^=0, (A6) 
where the subscript indicates that the function is evaluated at ?/ = 0. 



We have checked that (A4)-(A6) agree exactly with the covariant junction conditions, (20) 



and (22 ), specialized to this gauge. Further, they reduce to those of the standard DGP model 



in the strong coupling limit tt — 0, itlq — t- 0, in which the bulk scalar vr decouples. 
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FIG. 2: Example evolution histories in which no cosmological constant is present to drive cosmic 
acceleration. [Top panel] The deviation of the expansion history from that derived from stan- 
dard matter (for which 3H'^ / pm = 1). The blue and red curves each show consistent solutions 



to the modified Friedmann equation (32): one solution (red, thick line) recovers the standard ex- 
pansion history at early times and then undergoes accelerated expansion at late times; the other 
solution (blue, dotted line) has an expansion history entirely inconsistent with that of standard 
ACDM, with the vr field dominating the expansion at all eras, and undergoing heavily decel- 
erated expansion at late times. [Center panel] The evolution of the effective fractional energy 
density, = SttGptt/SH'^ , for the two solutions discussed above. For the accelerating solution, 
the phantom-like behavior in the matter era allows the vr field to dominate and drive cosmic ac- 
celeration at late times. The model is not physical, however, since as —7-2/3 one finds H ^ oo 
and a singularity occurs. [Bottom panel] A comparison of the effective equation of state for the 
expansion, Wes = — 1 — (2/3)dlnff/dlna, for the accelerating tt (red, full line) and fiducial ACDM 
(black, dashed line) scenarios. For the vr driven expansion histories, we use the numerical values 
Ho = 2.33 X 10""^ Mpc"^ tuq = 3.5 x 10"^^ Mpc^^ and ms = 4.4 x 10"3i]y[pc-i for which the 
maximum singularity occurs just after a = 1. 
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